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ABSTRACT 

We report the discovery of a z=2.76 radio galaxy, MG1019+0535, with an unusual spectrum. Lya, which 
usually dominates the spectra of high redshift radio galaxies, is very weak and the strongest lines in the spectrum 
are CIVA1549 and HeIIA1640. We speculate that dust is responsible for attenuating the Lya line. This object 
provides more evidence for dust formation at early epochs. 

The optical counterpart of the compact radio source has a double morphology which is orthogonal to the radio 
source axis and most likely due to the superposition of two distinct objects which may or may not be physically 
related. There is only one detected line from the brighter component (B) which is close in observed wavelength to 
the redshifted HeIIA1640 line from the z=2.76 galaxy (A). We discuss whether B is at the same redshift as A or 
whether it is a foreground galaxy at z=0.66. We also report the serendipitous discovery of a very red [R — K>7) 
field object. 

Subject heading: early universe — galaxies:active — galaxies: distances — galaxies:individual (MG1019+0535) 
— ISM:dust, extinction — radio continuum:galaxies 



1. Introduction 

The identification of the optical counterparts of faint radio 
sources remains one of the most effective methods of find- 
ing galaxies at high redshifts [e.g. Minkowski 1960, Lacy et 
al. 1994), primarily because all powerful radio galaxies have 
strong emission line spectra. The strongest line typically seen 
in the rest frame UV spectra of z>1.8 radio galaxies is Lyman 
a, a fact which led to early interpretations of high redshift ra- 
dio galaxies as protogalaxies. Although we now believe that 
the luminous Lya nebulae in high redshift radio galaxies are 
powered by their AGN rather than by stars, the discovery of 
strong Lya emitters at high redshift is still regarded as a pos- 
sible protogalactic signature (Partridge & Peebles 1967) and 
has prompted deep Lya imaging searches for primeval galax- 
ies. Thus far, these searches have all been unsuccessful (see 
Djorgovski et al. 1993 for a recent review). One possible rea- 
son for this is that dust suppresses the Lya emission from 
protogalaxies; indeed, even small dust optical depths can com- 
pletely extinguish the Lya emission (Hummer & Kunasz 1980, 
Neufeld 1990). 

Is there any evidence for dust at high redshift? Although 
there is weak statistical evidence for dust in damped Lya sys- 
tems {e.g. Fall & Pei 1993), more direct evidence has been 
found in the z=2.28 ultraluminous IRAS galaxy F10214+4724 
(Rowan-Robinson et al. 1991) and, more recently, in the 
z=2.34 steep spectrum radio galaxy TX0211-122 (van Ojik et 
al. 1994). Both these objects show very weak Lya lines and in 

^ Based on observations at the W.M. Keck Observatory 
Visiting Astronomer at Kitt Peak National Observatory, National Op- 
tical Astronomy Observatories, which is operated by Association of Uni- 
versities for Research in Astronomy, Inc., under contract by the National 
Science Foundation. 

Visiting Astronmer at the Very Large Array of the National Radio As- 
tronomy Observatory, which is operated by Associated Universities, Inc., 
under contract by the National Science Foundation. 



LI 



addition F10214 is a powerful infrared emitter. In this paper 
we report the discovery of a z=2.76 radio galaxy with weak 
Lya. It has unusual spectroscopic and morphological proper- 
ties and may be the most distant dusty galaxy discovered to 
date. 

Over the last few years we have been pursuing the opti- 
cal identification of a sample of 216 radio sources selected in 
a systematic manner from the SGHz MIT/Greenbank survey 
(MG I) of Bennett et al. (1986). (Several other optical studies 
of radio flux density selected samples are presently in progress; 
see McCarthy 1993 for a recent review.) Our sample is selected 
based on radio source angular size [d < 10") and radio spec- 
tral index (aj'^Qg^ ^ 0.75, F^ ~ z/~") which biases the sample 
toward flnding galaxies at high redshift. In addition, we only 
select sources which have simple [i.e. unresolved, double or 
triple) radio morphologies. The details of and results from 
our survey will be discussed in a future paper; some prelimi- 
nary results are described in Spinrad et al. (1993). The galaxy 
discussed here has the second highest redshift in our sample. 

Throughout this paper we use Ho = 50 km s~^Mpc~^ and 
qo = which implies an angular scale of 13.5 kpc/" at z=2.76. 
For Ho = 50 km s~^Mpc~^ and qo = 0.5, the corresponding 
scale would be 7.5 kpc/". 

2. Observations 

MG1019+0535 (J2000; 1016+0549 B1950) was observed on 
numerous occasions over the past 2 years at the KPNO May- 
all 4m and the Lick Observatory Shane 3m telescopes. The 
faint [R ~ 23) optical counterpart to this radio source was 
flrst identifled on images taken with the 4m in April 1990; 
deeper images were obtained in better seeing with the 4m in 
April 1992. Spectra taken with the 4m and RC spectrograph 
in November 1990 and with the lens/grism spectrograph on 
the Lick 3m in December 1990 and January 1991 showed the 
CIVA1549 and HeIIA1640 emission lines at a redshift z = 2.76, 
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but Lyman a A1216, normally by far the strongest UV spec- 
tral feature in high redshift radio galaxies, was only weakly 
detected. We considered this to be sufficiently unusual as to 
doubt the redshift until further observations obtained at Lick 
with the Kast spectrograph in February 1994 confirmed the 
redshift with a secure detection of CIII]A1909. These previous 
optical observations have largely been superceded by the Keck 
data presented below; we therefore discuss only these latter 
data here. 

We observed MG1019+0535 at the Cassegrain focus of the 
10m W. M. Keck Telescope using the Low Resolution Imag- 
ing Spectrometer (LRIS) on U.T. 1994 March 14 and 15. The 
detector is a Tektronix 2048^ CCD with 24/im pixels corre- 
sponding to a scale of 0"214 pix~^. Our observing run was 
one of the first Keck LRIS science runs and there were several 
problems with the instrument and detector: the readout noise 
of both amplifiers was excessively high (25e~ on the blue side 
and 65e~ on the red side) and the red side showed non-linearity 
at moderate count levels which resulted in difficulties in flux 
calibration at wavelengths longer than 7300A. The typical see- 
ing during our observing run was 0"8 - l"0. 

Three direct images of MG1019+0535 (total i<.xp = 1300s) 
were obtained through a Cousins R (3mm OG570 + 3mm KG3; 
Xeff fs 6500A, AA fs 1400A) fllter. The images were corrected 
for overscan bias, flat flelded using a median sky flat and coad- 
ded. The effective seeing in the coadded frame is l"0 and the 
3(7 detection limit in a 3" diameter aperture is R f«27.2 mag. 
The imaging data were calibrated using observations of the 
Landolt (1992) standard star fleld PG1525-071 and the spec- 
trophotometric standard G193-74 (Oke 1990). 

Five spectra (total iexp = 7500s) of MG1019+0535 in the 
wavelength range 3900— 88OOA were obtained through a 1" slit 
in PA=36° using a 300 1/mm grating (A^ = 5000A, resolution 
FWHMfsloA). The slit was aligned on the object by using 
a blind offset from a bright offset star which lies 26"6 west 
and 28"0 north of the MG source. The slit PA was chosen to 
trace the axis through the two components of the galaxy; it 
was also parallel to the average parallactic angle during our 
observation. Hence, although the spectroscopic data were ob- 
tained under non-photometric conditions, we believe our rela- 
tive spectrophotometry to be accurate for A < 7300A. The data 
were corrected for overscan bias and flat flelded using internal 
lamps taken immediately following the observations. Spec- 
tra of the two components were extracted using parallel linear 
models for the centroid variation along the dispersion axis and 
were flux calibrated using observations of the spectrophoto- 
metric standards G193-74 (Oke 1990) and IIZ44 (Massey et 
al. 1988). The standard stars were observed both with and 
without an OG570 fllter in order to correct for the second or- 
der light contaminating the wavelength region A > 7500A. 

The fleld was imaged in the Ks band (1.99-2.32/im) us- 
ing the Infrared Imager (IRIM) on the KPNO 4m on UT 
1994 February 25. IRIM employs a 256x256 pixel HgCdTe 
NICMOS-3 array with a scale of 0!'603/pixel at the 4m F/15 
Cassegrain focus. Thirty frames were obtained, each an av- 
erage of three 20 second exposures, with the telescope off- 
set between frames in a non-redundant pattern. Array re- 
sponse variations were corrected using dome flat exposures. 
Sky frames were created for and subtracted from each image 
using a object-masked combination of eight temporally adja- 
cent exposures. The resulting frames were spatially oversam- 
pled by a factor of four, co-registered to the nearest sub-pixel. 



and summed. The total exposure time in the resulting im- 
age is 1800 seconds, and the image PSF has a FWIIM=1"1. 
Sky conditions were judged to be photometric based on min- 
imal frame-to-frame zeropoint variations and the excellent 
quality of the sky subtraction; the zeropoint uncertainty is 
small compared to the random error in the photometry on this 
faint source. Observations of faint UKIRT infrared standards 
(Casali & Hawarden 1992) provided photometric calibration, 
under the assumption that K=Ks for these stars (reasonable 
given their nearly uniformly zero H—K colors). The reduction 
of all our imaging and spectroscopic data were done using the 
IRAF package. 

MG1019+0535 was observed in snapshot mode using the 
VLA in the A-array conflguration on 12 August 1991 at 
8.44GIIz and 1.4GIIz. The exposure time in each band was 
roughly 7 min. We observed the VLA calibrator 0922+005 for 
phase calibration and 3C286 for flux calibration. Calibration 
and analysis of the radio data were carried out using the AIPS 
package. The data were self calibrated in three iterations, the 
flrst two iterations using only phase self calibration and the 
last using amplitude self calibration. 

3. Results 

Figure 1 shows a detail from our Keck i?-band image of 
MG1019+0535. The galaxy has two components separated 
by l"5 in PA=36° of which the southern, fainter component 
('A') is probably associated with the radio source. In §4.2, we 
discuss the interpretation of this double morphology in some 
detail. Here, we simply note that the i?-band image is well 
flt by the superposition of two distinct model galaxies with 
elliptical isophotes. A measurement of the R magnitude of 
each component was obtained after subtracting out the best 
flt model to the other component. We present the photometry 
in Table 1, including a pure-continuum R magnitude for ob- 
ject A after removing a 17% emission line contribution derived 
from our spectroscopy. 

A panel showing both the R and Ks images is shown in 
Figure 2 (Plate). The K image of MG1019+0535 is crudely 
consistent with the two component morphology of the R image. 
Although the signal-to-noise ratio in the K frame is too poor 
to allow accurate division of the total flux between the two 
components A and B, we have attempted a rough estimate by 
using the two elliptical model constructed for the R image. We 
scaled the flts to the two i?-band components independently 
and subtracted them from the K image until the residuals 
were minimized. The resulting magnitudes for the individual 
components are listed in Table I with an estimate of the errors. 
The Ks band is free of strong emission lines at the redshift of 
the radio galaxy. 

The VLA map of the radio source is shown in Figure 3. 
The radio source has a triple structure with two bright lobes 
(/I.'41ghz = 24mJy,/8--*(3Hz = 27.5mJy) separated by about 
l"3 and a faint core. The spectral index integrated over the 
radio source is af;|^*^^^ = -1.0. The radio axis PA of 103° is 
roughly orthogonal to the line joining the two optical compo- 
nents. This is unusual for high redshift radio galaxies which 
generally show alignments between their rest frame UV mor- 
phologies and their radio axes (McCarthy et al. 1987, Cham- 
bers et al. 1987). Table 1 lists the radio and optical positions 
and flux densities of the various components of MG1019+0535. 

The cross in Figure 1 marks our best estimate of the opti- 
cal position of the radio core based on 5 astrometric stars in 
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the field. Using different combinations of astrometric stars re- 
sulted in slightly different positions for the core position (with 
slightly larger errors) but these all appear to scatter close to 
the position of the optical component A. Therefore we con- 
clude that the radio source is most likely to be associated with 
component A. If the core coincides with the optical compo- 
nent A, then the extended radio structure is entirely confined 
within the optical extent of the galaxy. Component A may 
be slightly extended in the east-west direction implying that 
it may be aligned with the radio source axis. 

The LRIS spectra of the two optical components are shown 
in Figure 4. The central wavelengths, fluxes, equivalent widths 
and full widths at half maximum of each emission line in the 
observed frame are listed in Table 2. The fluxes listed are those 
measured in a l"xl"7 aperture. Component A has strong and 
narrow associated line emission from Lya, NV, CIV, Hell and 
cm] at a (weighted) mean redshift z=2.765±0.002 (Fig. 4a). 
The Lya line is redshifted with respect to the CIV, Hell and 
CHI] lines by ~500kms~^; this is similar to offsets observed 
in other high redshift radio galaxies and quasars. In addition, 
there is a tentative detection of OVIA1034 at Ao6s = 3903 A. 
This feature is at the very edge of our spectral window where 
the instrumental sensitivity is decreasing and therefore should 
be regarded as uncertain. 

The CIV, Hell and CHI] lines are broader than the Lya and 
NV lines. For reference, the night sky lines of [OI]A5577 and 
[OI]A6300 have a FWHM fsloA in our spectrum. This implies 
that the Lya and NV lines are unresolved, whereas the Hell 
line is broadened by a velocity dispersion of ~400kms~^. 

Component B has a single, unresolved, moderately strong 
line at 6192A, separated by only 17A in observed wavelength 
from the redshifted Hell line in component A. Its continuum 
spectrum appears to have a discontinuity at A fs 4500A below 
which there is little detectable flux. The spectrum appears to 
rise longward of A « 6300A. The drop in the flux at A > 7600A 
is not real and probably due to instrumental effects [e.g. prob- 
lems with the red side CCD readout amplifler, the decreasing 
CCD sensitivity at long wavelengths and the increased read 
noise on the red side). The optical spectrum of B appears red- 
der than that of A, although its i? — if color is slightly bluer. 

Since only one strong line is detected from component B, 
we cannot be certain of its redshift. Three possibilities are 
(1) that the line is Lya, (2) the line is [OII]A3727 which is 
one of the strongest lines in low redshift star forming galaxies, 
and (3) it is HeIIA1640. The last possibility is based primar- 
ily on the proximity of the line in B to the HeIIA1640 line 
in A. In addition, whether or not B is physically associated 
with component A, the unusual spectrum of A still requires an 
explanation. These issues are discussed in §4.2. 

4. Discussion 

4.1 The Unusual Spectrum of Component A 
AND THE Evidence for Dust 

The most unusual aspect of the spectrum of component A 
is that the Lya line is very weak compared to that observed 
in most high redshift radio galaxies. Table 3 compares the 
various line strengths (relative to HeIIA1640) with the corre- 
sponding ratios in a few 'canonical' objects. The composite 
radio galaxy spectrum was constructed from 3CR and 1 Jy- 
class galaxies with 0.1<z<3 by McCarthy (1993). The UV 
spectrum of NGC1068 was obtained by Kriss et al. (1992) us- 
ing the Hopkins Ultraviolet Telescope. The only other reported 



examples of high redshift objects with suppressed Lya are the 
dusty, ultraluminous IRAS galaxy F10214+4724 at ^=2.286 
(Rowan-Robinson et al. 1991, Elston et al. 1994), and a z=2.34 
radio galaxy TX0211-122 (van Ojik et al. 1994). TX0211-122 
has been argued to be a dusty galaxy based on the similarity 
of its rest frame UV spectrum to that of F10214+4724. For 
comparison. Table 3 also lists the line ratios predicted by the 
photoionization code CLOUDY (Ferland 1994) for two values 
of the ionization parameter [U = Q[H.°)/4irr'^Ne) of —1.8 and 
—2.8. In both cases a power law input spectrum of S^, ~ z/~^-® 
and a density of 10 cm~^ was used. 

Although the CIV/Hell/CIH] ratios observed in object A are 
similar to those in McCarthy's composite radio galaxy spec- 
trum, the NV/Hell and Lya/Hell ratios in A are respectively 
factors of 2 and 10 times lower than those in the composite 
spectrum. The Lya/Hell ratio varies only slowly with ioniza- 
tion parameter [e.g. Ferland & Netzer 1983) and therefore it 
is unlikely that the observed ratio is due to differences in U, 
although varying the shape of the input spectrum may change 
the observed ratio. It is interesting to note that the fact that 
the CIV/Hell/CIH] ratios are all roughly unity implies that 
the ionization parameter is about 0.01 [e.g. Ferland & Netzer 
1983), which is similar to that derived for most radio galaxies. 
We do not see the anomalously high NV strength observed in 
the spectra of F10214+4724 and TX 0211-122 (van Ojik et 
al. 1994). 

In addition to component A's unusual line ratios, its Lya lu- 
minosity is also intrinsically weak (i/Lya = 1-6 x 10'*^ergs~^). 
McCarthy (1993) presents an empirical correlation between 
emission line luminosity and 1.4 GHz radio power. Using 
canonical values for the Lya/ [OH] ratio of 5 to 8 (also from 
McCarthy 1993) to convert his relationship for our purposes, 
and assuming negligible contribution from HeIIA1215 to the 
observed line flux from MG 1019+0535 (see below), we flnd 
that its Lya luminosity is a factor of 30 to 50 below the em- 
pirical correlation. We note, however, that the scatter in the 
i[oii] vs. Pi.4GHz relation is very large, being almost an or- 
der of magnitude. The Lya line in MG 1019+0535 does not 
appear to be spatially extended in our two-dimensional spec- 
tra, either in our Lick/KPNO data with the slit oriented along 
the radio axis (approximately PA=90°) or in our Keck spectra 
with the slit oriented along the line joining the two compo- 
nents (PA=36°). This is in sharp contrast to most powerful, 
high redshift radio galaxies, where the Lya emission is both 
strong and spatially extended (generally along the radio axis) 
over many tens of kpc [e.g. 3C294 McCarthy et al. 1990). 

The most straightforward method of destroying Lya is us- 
ing dust. Since Lya is a resonance line, multiple scatterings 
by ground-state hydrogen atoms signiflcantly increase the path 
length traversed by a Lya photon in a nebula thus increasing 
the probability of absorption of the photon by dust (Hummer 
& Kunasz 1980, Neufeld 1990). In MG1019+0535, a modest 
E[B-V) of 0.3 or 0.4 and a Galactic extinction curve will cor- 
rect the observed NVA1240/HeIIA1640 ratio to 0.44 - 0.51, a 
value similar to that measured for the 'average' radio galaxy 
(Table 3). This E[B-V) estimate assumes that the dust ex- 
tinction is mostly external to the line emitting region [e.g. a 
foreground dusty screen), but if the relative velocity between 
this dusty gas cloud and the line emitting gas cloud is negligible 
and the velocity gradients within the two clouds are compa- 
rable, then most of the emitted Lya photons will be trapped 
in the dust cloud. A similar quantity of dust intermixed with 
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the line emitting gas could also be responsible for attenuat- 
ing the Lya line as long as the velocity gradients in the line 
emitting gas are small so as to prevent most of the Lya pho- 
tons from escaping. This amount of extinction is more than 
sufficient to completely suppress the Lya emission relative to 
IIeIIA1640. For example, in the recent calculation by Chen & 
Neufeld (1994), it is found that for even small total absorption 
optical depths, the fraction of photons that escapes in the cen- 
ter of the Lya line can be more than an order of magnitude 
lower than the fraction that escapes at, for example, the po- 
sition of the NV line 24A away from the Lya line center. In 
fact, it is extraordinary that we observe Lya at all in distant 
galaxies. The geometry of the line emitting region and the spa- 
tial distribution of dust and dust to gas ratio in most objects 
must conspire to permit the ionizing radiation to reach a rela- 
tively dustless environment. Hence, although it is possible that 
there are other more involved methods of suppressing the Lya 
line, the simplest method which also is consistent with the ob- 
served spectrum is extinction by dust. We therefore conclude 
that MG1019+0535 is probably a dusty galaxy. 

Although the determination of the extinction using the 
NV/Hell and NV/CIV ratios is uncertain due to the intrin- 
sic variation of these ratios that is observed in radio galaxies, 
the estimate of E(B— V)=QM is comparable to other reddening 
estimates for high redshift radio galaxies. For example, Mc- 
Carthy et al. (1992) obtained an estimate of E{B-V)piQ.l for 
two z=2.4 radio galaxies by using the Lya/Ha line ratio. Both 
of these galaxies have very large Lya/CIV ratios (>5) and it is 
therefore likely that the reddening is larger in MG1019+0535 
which has Lya/CIV fs 1. 

We can estimate an upper limit to the continuum reddening 
by requiring the intrinsic continuum to be no bluer than flat 
spectrum (i.e. /^ ~ z/'^,/c < 0), the bluest color observed for 
distant radio galaxies thus far (e.g. Eisenhardt and Dickinson 
1992, Bales et al. 1993). With this assumption, the maximum 
allowable extinction would be E(B— Vj<0.43. If the contin- 
uum and lines have a common extinction, the reddening is con- 
strained to be 0.3 < E[B—V) < 0.43. However, the narrow line 
emission and continuum components may originate in physi- 
cally distinct regions as is probably the case for F10214+4724 
where the extinction derived from the emission lines would 
deredden the continuum to be bluer than any observed galaxy 
or AGN spectrum (Elston et al. 1994). 

If we assume that the extinction in MG1019+0535 is 
E(B— Vj=0.3, then for the mean Galactic extinction curve 
of Savage & Mathis (1979) the extinction coefficients for 
the R and K fllters in the rest frame of MG1019+0535 are 
Aij,^=2.766 = 8.08^(5 - V) and ^^^,^=2.766 = 2.92^(5 - V). 
The dereddened continuum colour (after correcting for line 
contamination to the R photometry) would then be [R — 
K)=2.i&. 

An observation of the HeIIA4686 line (which is redshifted 
into the H band at 1.5454/im) would provide an indepen- 
dent measure of extinction. The Case B recombination values 
for this ratio lie between 5.4 and 9.0 for 3000K<T<50000K 
and 100cm-2< re,. < lO^'cm-^ (Hummer & Storey 1987). 
Hence the predicted flux of the HeIIA4686 line for the extreme 
(T'=50000K, ree = 10''cm-2) case is « 10-^'' erg cm-^ s-^; this 
is a lower limit if dust extinction is important. 

Other evidence for dust in high redshift galaxies has been 
presented by McCarthy et al. (1992) and Bales & Raw lings 
(1993) who measured the Lya/Ha and Lya/H/3 ratios in a 



few high redshift radio galaxies. These line ratios are mea- 
sured to be much lower than predicted by photoionization 
models, implying the existence of dust. The one high redshift 
object which shows convincing evidence for dust is the IRAS 
source F10214+4724; this object has suppressed Lya, a large 
far infrared luminosity presumably due to dust emission, and 
large amounts of molecular gas (Rowan-Robinson et al. 1991, 
Solomon et al. 1992). If MG1019+0535 were to have the same 
infrared luminosity as F10214+4724, its larger redshift would 
cause the measured infrared flux densities to be ~ 1.5 - 2 times 
fainter, rendering it undetectable in the IRAS database. 

Table 3 shows that the HeIIA1640 emission line is compa- 
rable in strength to Lya (Fig. 4a, Table 2). The HeHA1640 
line is the "Balmer" a (3— >2) transition of the Hell ion. The 
corresponding Hell Balmer j3 transition is at 1215.17A and 
therefore contributes to the measured hydrogen Lya flux. In 
most known high redshift galaxies the HeIIA1640 line is very 
weak relative to Lya and this contribution is neglible. How- 
ever, in the present case the observed HeIIA1640 line is very 
strong and the contribution needs to be considered. For Case 
B recombination, the HeIIA1215/HeIIA1640 ratio ranges be- 
tween 0.27 and 0.33 for temperatures 5000K<T'<20000K and 
densities lOOcm"^ < Ne < 10'*cm~^ (Hummer & Storey 1987, 
Osterbrock 1989). As a result, the "true" Lya/HeIIA1640 ratio 
in MG1019+0535 probably lies between 0.67 and 0.72. 

Note that the strong HeIIA1640 line implies that the nebular 
gas has a high degree of ionization which lends some credence 
to our tentative identiflcation of the OVI line. If real, the 
OVI/Hell ratio of ~ 0.4 is similar to the ratio of 0.6 measured 
in the composite radio galaxy spectrum of McCarthy (1993). 
In comparison, the HUT spectra of NGC4151 and NGC1068 
show OVI stronger than HeIIA1640 (Kriss et al. 1992a, 6). 

4.2 Morphology and Photometry 

The most striking aspect of Figure 1 is the double morphol- 
ogy of the system which is reminiscent of low redshift double 
galaxy systems (e.gr. 3C89, 3C169.1, 3C433). However, the fact 
that we detect only one strong line in the spectrum of compo- 
nent B raises the question as to whether (i) A and B are at the 
same redshift (z=2.76) and separated spatially by f«20 kpc (in 
projected distance) and in velocity by only f«800kms~^ or (u) 
whether the two components are unrelated and at two different 
redshifts. In this latter case, the strong emission line seen in 
the spectrum of component B is most likely to be [OII]A3727 
implying that B is at z=0.66. 

The K magnitude for component A alone is faint for its red- 
shift, about 1 magnitude less luminous than the median value 
implied by the K vs. z diagrams in McCarthy 1993 and Bales 
et al. 1993. Given the observed scatter in radio galaxy K mag- 
nitudes at these large redshifts, however, this deviation from 
the typical values is perhaps not unusual. If the radio galaxy 
were instead considered to consist of components A plus B 
together, its if-luminosity would be typical for its redshift. 
The color [R — K m 3.9 for component A alone, R — K m 3.1 
for the entire system A+B) is not unusual for radio galax- 
ies at z>2. If we apply the extinction correction implied by 
E[B — V)= 0.3 to the continuum of object A, its K luminos- 
ity would be consistent with the K vs. z prediction, although 
the resulting corrected R magnitude would then be ~1 to 2 
magnitudes brighter than the norm. 

In what follows, we discuss the evidence for the 'discrepant 
redshift' and 'same redshift' cases. 
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4-8.1 A Chance Superposition.... 

The main argument for different redshifts for A and B is that 
while component A shows emission lines from Lya, CIV, Hell 
and CIII], component B has only one line. In our experience, 
the detection of a single strong line over a large searched range 
in wavelength usually implies that the line is either Lya or 
[OII]A3727. It is unlikely that the single line from B is Lya 
because (i) there is no drop in the continuum shortward of the 
line which would be expected due to Lya forest absorbers and 
(u) there is some detectable continuum shortward of 4650A 
which would correspond to 912A in the rest frame. On the 
other hand, it is entirely plausible that the line is [OII]A3727 
at z=0.66. In this case, the luminosity of the line is large 
(ijoiii « 8 X 10"*° ergs s~^), but not unprecedented for star 
forming galaxies. The rest frame equivalent width of the line 
is VKjoii] Si 26A which is well within the range of equivalent 
widths measured in the integrated spectra of nearby galaxies 
{e.g. Kennicutt 1992a, Colless et al. 1990, 1994). In addition, 
the single emission line from B is unresolved, which implies 
that the width of the line is probably less than 240 kms~^. 
Hence it is unlikely that B is a luminous AGN. 

If B is at z = 0.66, then the rising continuum redward of 
Xrest f«3800A [Xobs f«6300A, see §3 and Fig. Ah) is typical of 
star forming galaxies, in which it is dominated by the con- 
tinuum from intermediate age stars (in contrast to the "true" 
4000 A break feature seen in old stellar populations). The break 
amplitude for B is Dy(4000)=1.56, quite typical for star form- 
ing galaxies at z=0.6 (Bruzual & Chariot 1993). The R — K 
color of object B is also consistent with that of an actively star 
forming galaxy (e.g. a late-type spiral) at z = 0.66. 

Ignoring the absence of other strong emission lines for the 
present (see §4.2.2), it is of note that almost all of the galax- 
ies in the Kennicutt (1992a, 6) galaxian spectral atlas with 
W'^loii] ^ 25 are very late-type galaxies: except for 1 SO which 
is a Seyfert, they all either have Hubble types of Sc or later, 
or are peculiar or interacting galaxies. If B is at z=0.66, then 
the i?-band filter samples the continuum light between f«3300- 
4800A. The resolved morphologies of most late-type galaxies 
are more diffuse and much less nucleated at these wavelengths 
than the morphologies of early-type galaxies [e.g. Schweizer 
1976). Although our ability to distinguish small scale mor- 
phological details at this redshift is degraded by our resolution 
(the diameter of our seeing disk corresponds to f«7.5 kpc at 
z=0.66), simulations show that the overall surface brightness 
profile of component B is better fit by an exponential disk con- 
volved with the seeing than a deVaucouleurs profile. Therefore 
the morphology of component B is consistent with that of a 
z=0.66 late-type spiral. 

There are two arguments against the low redshift hypothe- 
sis for B. Firstly, we do not detect other strong lines (such as 
[OIII] and H/3) that are commonly observed in the spectra of 
low redshift emission line galaxies. This is discussed in further 
detail in §4.2.2. Secondly, low redshift emission line galaxies 
are not known to have a spectral discontinuity at Xrest f«2700A 
{Xobs f«4500A; see Fig. 4b). The UV spectra of early-type 
galaxies can be steeply declining in this region. However, the 
optical spectrum, the optical morphology and the optical-IR 
color suggest instead that B would be a late-type spiral which 
would therefore have a much flatter UV spectral energy distri- 
bution. 

If the double morphology of the system is due to a chance 
superposition, is this an extremely unlikely event? Let us con- 



sider the probability that we would flnd a chance superposition 
between a radio source and an unrelated fleld galaxy along the 
same line of sight. There are 216 radio sources in our sample, 
of which we have imaged 100 to date. The surface density of 
galaxies with R magnitudes between 18.0 and 23.5 is approx- 
imately E = 3xl0~^ arcsec"^ (Metcalfe et al. 1991; galaxies 
brighter than R=1S contribute insigniflcantly to the counts). 
Therefore, the a priori expected number of superpositions be- 
tween unrelated fleld galaxies and radio sources within d" of 
one of the radio sources in our sample is crudely IOOEtt^^, or 
~2 for d=l'.'5. Hence it is not improbable that we should flnd 
such a superposition. Here we have neglected the angular cor- 
relation function of galaxies under the assumption that chance 
projections should be uncorrelated; including it into the calcu- 
lation roughly doubles the probability. Gravitational lensing 
by the foreground galaxy could enhance the likelihood for such 
a projection, but we see no obvious evidence for this in the case 
of MG 1019+0535. In fact component A is rather underlumi- 
nous at K compared to other radio galaxies at similar redshift. 

Using the redshift distribution of faint fleld galaxies from 
combined / < 22.5 samples of Lilly (1992) and Tresse et 
al. (1993) (which should be roughly equivalent to that for 
galaxies at i? < 23), the probability that a given fleld galaxy 
would have a redshift placing [OII]A3727 within 17A of a 
strong line in the radio galaxy (Lya, CIV, Hell or CHI]) is 
0.036. Of the galaxies in those surveys for which the spec- 
troscopic observations cover the [OH] line, ~78% have de- 
tected emission. Thus the combined expectation value for 
the observed superposition using the joint angular proximity, 
wavelength coincidence, and emission line fraction criteria is 
roughly 2 x 0.036 x 0.78 = 0.056. 

In summary, it is plausible that B is a foreground object, 
although it would have somewhat unusual spectral properties. 
The similarity of the emission line wavelengths would, in this 
scenario, be a somewhat remarkable coincidence, but not an 
inflnitesimally likely one. 

4.S.S ....Or A Double Galaxy System? 

An important argument against the z = 0.66 scenario for 
component B is the absence of [OIII] A5007, 4959 and H/3 emis- 
sion, which should appear in our spectral window. If B is at 
z=0.66, the rest frame equivalent width of the strong emission 
line in B is 26A. Although this is fairly typical of widths ob- 
served in low redshift galaxies [e.g. Kennicutt 1992a, Colless 
et al. 1990,1994), galaxies that exhibit similar VKroii] are usu- 
ally star forming galaxies or AGN. In either case the [OIII] 
lines would also be strong and easily detectable. For star 
forming galaxies, the observed [OIII] /[OH] ratio exhibits a 
large range from 0.1 to 1 and in particular, for galaxies which 
have VKjoii] > 25A, the ratio is usually observed to be greater 
than 0.2 {e.g. Kennicutt 1992a). For most AGN other than 
LINERs, this ratio is typically greater than 1 [e.g. Baldwin, 
Phillips & Terlevich 1981). However, based on the spectrum of 
B, the limit on [OIII]/ [OH] <0. 03 implies that B is an unusual 
object if it is at the lower redshift. 

The spectrum of a z=0.544 fleld object that we ob- 
tained serendipitously with the spectrum of MG1019 pro- 
vides a useful comparison. This serendipitous object 
lies about 44" east and 60!'5 north of MG1019+0535. 
Its spectrum shows strong [OII]A3727 emission [fobs = 
4.3 X 10~^'^erg s~^ cin~'^ ,Wx,rest = 65 A) but also exhibits 
strong [OIII] A4959, 5007 and H/3 emission with the ratios 
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[OIII]A5007/[OII] = 1.33 and H/3/[0II] = 0.47. Although the 
serendipitous galaxy was not centered on the slit and therefore 
the line ratios may only be representative, it is important to 
note that we were easily able to detect the [OIII] and II/3 lines 
in this galaxy but not in component B. 

An additional clue is obtained from comparing the flux ra- 
tio across the Lya line in the two components. For compo- 
nent A, the Lya continuum (forest) discontinuity (/x(1300 — 
1400)//a(1100- 1200)) is roughly 1.5±0.1. However, B has 
very little flux below Xohs ~ 4500 (Fig. 4b). If B is at z=2.76, 
this translates into a 'Lya forest' continuum discontinuity of 
about 4.5±0.5. This is extremely large; typical values of the 
forest discontinuity at z=2.7 lie between 1.2 and 1.7 (Schnei- 
der, Schmidt & Gunn 1991). Hence, if this break is real, it 
can only be explained by invoking some strongly absorbing 
gas such as would result in a strong damped Lya system or a 
BAL-type system. Unfortunately we can neither conflrm nor 
rule out whether this break is an instrumental effect; a com- 
parison with the serendipitous z=0.544 galaxy is inconclusive 
since its spectrum has very little detected continuum below 
5000A. However, as discussed in §4.2.1, a break of this ampli- 
tude is not usually seen in emission line galaxies. 

If both A and B are at z=2.76, then the strongest UV line 
in the spectrum of B is HeIIA1640 which is unusual (perhaps 
unprecedented). However, there is no a priori physical reason 
why this situation could not occur. Simulations with CLOUDY 
suggest that the CIV/Hell ratio is fairly sensitive to the ioniza- 
tion parameter. Speciflcally, decreasing the ionization param- 
eter from \og{U) = -1.8 to \og{U) = -2.8 for an /^ ~ i^-^-^ 
input spectrum on solar abundance gas changes the CIV/Hell 
ratio from 1.13 to 0.01 (Table 3). If we assume that the sole 
source of ionization in B is due to the AGN nucleus in A, then 
B would experience a highly diluted radiation fleld relative to 
that experienced by the ambient medium in A. Alternatively, 
the ionization parameter may be decreased in B relative to A 
by invoking a higher particle density for the ionized gas in B 
relative to A. The Lya line from B is easily suppressed by dust. 
Although suppressing the CIV relative to Hell requires lower- 
ing the ionization parameter an order of magnitude from the 
'canonical' model, U must be decreased an additional order of 
magnitude (log({7)<-4) to suppress the CHI] line relative to 
HeH. 

Given that the Lya line implies the presence of dust, it is 
suggestive that the double optical morphology of the system 
resembles that of nearby dusty radio galaxies such as Centau- 
rus A (Baade & Minkowski 19546, Burbidge & Burbidge 1959) 
and Cygnus A [e.g. Baade & Minkowski 1954a). Could it be 
that we are observing a dust lane in this high redshift galaxy? 
There are two arguments against this suggestion. Firstly, as 
noted above, we can model the observed structure well by us- 
ing two separate objects with undistorted elliptical isophotes. 
The morphology is therefore consistent with the superposition 
of two separate, undistorted galaxies (without specifying their 
relative redshifts). This implies (invoking Occam's razor) that 
a dust lane is unnecessary in order to explain the morphology. 
Secondly, the optical astrometry places the radio core position 
very near component A. If the double optical system were in- 
deed a single galaxy with a dust lane, one might naively expect 
the radio source to be centered between the two components 
and orthogonal to the radio axis, as is observed in Centaurus 
A (Wade et al. 1971, Meier et al. 1989). 

If the two components are instead separate objects at the 



same redshift, then their projected separation is only 20 kpc 
and their observed difference in velocity is 800 kms~^. This is 
of the order of cluster velocity dispersions. In fact, even larger 
relative velocities have been measured for double galaxies that 
are seen in association with radio sources. One such exam- 
ple is 3C 169.1 (z=0.633), in which the two components of the 
double galaxy system are separated by only 3" (27 kpc pro- 
jected distance) and have a relative velocity of 1400 kms~^! 
The observed morphology and kinematics of low redshift radio 
galaxies with double components have been argued by vari- 
ous authors to support the hypothesis that these systems are 
merging galaxies and that the AGN activity has been trig- 
gered by the merger {e.g. Heckman et al. 1986). With such 
a small projected separation and (if they are at similar red- 
shifts) moderate relative velocity, it is possible that the the 
two optical components of MG1019+0535 are in the process of 



mergmg. 



4.2.3 An Extremely Red Object 



A visual comparison of the R and K images of the 
MG1019+0535 fleld reveals an extremely red object (D in 
Fig. 2), which is fs 10" W of the radio source. While it is easily 
detected in the infrared image (ii'=18.26±0.12 in a 4" diame- 
ter aperture), it is undetected in our deep Keck R image. The 
object is nonstellar: it is spatially extended in K with a decon- 
volved FWHM fs 0"5. Although its proximity to the brighter 
galaxy C makes accurate photometry difiicult in the R image, 
a conservative estimate of its color is R — K > 7. Without 
invoking dust extinction, such red colors are most naturally 
explained as resulting from the large K-corrections expected 
for relatively unevolved, early-type galaxies at z > 1. The 
presence of a presumably foreground galaxy (C) close to the 
line of sight to D, however, could result in additional redden- 
ing. The bright K magnitude of the object makes it unlikely 
that it is at the redshift of MG 1019+0535. 

Similar very red objects have been discovered close to sev- 
eral powerful AGN at large redshifts. For example, McCarthy, 
Persson & West (1992) reported the discovery of very red ob- 
jects [R — K > 6, if ~ 18) near the radio galaxies 0406-244 
(z=2.427) and 2025-218 (z=2.630). The fleld of B2 0902+34 
(z=3.395, Lilly 1988) contains a, K = 18.3, R-K = 6.1 object 
which lies 15" from the radio galaxy (Eisenhardt & Dickin- 
son 1992). Graham et al. (1993) discovered two red objects 
with identical magnitudes and colors [K = 18.4, R — K = 
6.2) that lie 10" and 20" away from the high redshift galaxy 
4C41.17 (z=3.80. Chambers, Miley & van Breugel 1990). The 
two I - K K. 6.5 [K = 18.4 & 18.7) objects reported by 
Hu & Ridgway (1994) are in the fleld of the z=3.79 quasar 
PC 1643+4631A; one of the objects is less than 20" from the 
quasar. 

Infrared surveys of blank flelds have also discovered very red 
galaxies with bright K apparent magnitudes (cf. Elston, Rieke 
& Rieke 1991). The fleld surface density of the reddest objects 
is <0.01/n' based on a recent reanalysis of the Hawaii Deep 
if-band Survey data by Hu & Ridgway (1993). In compari- 
son, the detections summarized above imply that the surface 
density of very red objects in the flelds around high redshift 
radio galaxies is unusually large: > 1/n'. Obviously a more 
thorough statistical study of the presence of very red objects 
in the flelds of high redshift AGN with a proper treatment of 
the fleld population is required before drawing any conclusions 
on the basis of the few flelds mentioned here. If these objects 
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are z>l ellipticals, then an excess of these objects in the fields 
of high redshift AGN may imply that gravitational lensing is 
biasing our samples of the most distant radio galaxies. On 
the other hand, if they are at the redshift of the AGN and 
their emission is primarily starlight, these objects would pro- 
vide stringent constraints on the epoch of galaxy formation. 
The bright K magnitudes of the red objects, however, imply 
that this latter possibility is unlikely. 

5. Conclusions 

We have discovered a dusty z=2.76 radio galaxy as part of 
our ongoing study of radio sources from the MIT-Greenbank 
radio survey. The optical image shows two components sepa- 
rated by 1"5. Although these may be physically related, there 
is circumstantial evidence that one of the components is a 
foreground object at z=0.66. Whether or not the two com- 
ponents are related, the spectrum of the high redshift object 
is extremely unusual: it has very weak Lya and an observed 
Hell/Lya ratio of fsl. The most straightforward method of 
suppressing Lya relative to Hell is by dust extinction of the 
resonance line; therefore MG1019 is a good candidate for a high 
redshift dusty radio galaxy. The strength of the HeIIA1640 
emission line implies that HeIIA1215 might be an important 
contribution to the observed H Lya line strength. 

We have discussed the likelihood that the two optical com- 
ponents are at the same redshift or at two different redshifts. 
If they are at different redshifts, then the spectrum of B is un- 
usual with a moderately luminous [OII]A3727 line and a small 
[OIII]/[OII] ratio. However, the gross features of the spectrum 
[i.e. the strength of the single line and the shape of the con- 
tinuum) appear to weigh in favor of the conclusion that com- 
ponent B is at z=0.66. If, however, they are both at z=2.76, 
then component B has an even more unusual spectrum with 
HeIIA1640 being the only strong UV line. We speculate that it 
is possible for such a situation to arise if the line emitting gas 
in B is ionized by the diluted radiation field of the AGN in A. 
In addition, the observed morphology may imply that we are 
observing an interacting system at z=2.76. 



The discovery of dusty systems (or even galaxies with weak 
Lya) at high redshift has several important implications. 
Firstly, it implies that population synthesis modelling of the 
broad band colours of high redshift radio galaxies may not be 
straightforward since these models so far do not account for 
any reddening of the continuum by dust. Secondly, it is im- 
portant to note that there are only three known high redshift 
(z>1.8) radio galaxies with weak Lya. This implies that ei- 
ther dusty high redshift radio galaxies are very rare objects and 
most high redshift radio galaxies are largely dust free, or the 
dusty objects are missing from our high redshift samples due to 
selection and observational biases. Thirdly, it provides further 
evidence that dust extinction may be the reason that primeval 
galaxy searches designed to look for strong Lya line emission 
have been largely unsuccessful. Finally, the existence of dust 
at high redshifts raises the question of how and when this dust 
was formed. The presence of dust implies an early episode of 
star formation. This suggests that it is more appropriate to 
interpret the spectral energy distributions of high redshift ra- 
dio galaxies in terms of an underlying old population and a 
superimposed burst population (cf. McCarthy 1993) than in 
terms of a single burst of star formation. 
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Figure Captions: 

Figure 1: Keck LRIS i?-band image of MG1019+0535. The cross marks our best estimate for the position of the radio core. The 
scale bar of 1" corresponds to 13.5 kpc at z=2.765 and 9.3 kpc at z=0.66 (IIo=50kms~^ Mpc~^, qo=0). The position of component 
A is roughly aiggo = 10''16'"56!'8, ^iggo = +5°49'39". The R magnitudes of A and B are 23.7 and 23.0 respectively. 

Plate 2: Co-registered R and Ks images of the MG1019+0535 field. The extremely red object marked D in the Ks image is 
undetected in our deep Keck R image [R — K >7) where its location is denoted by the cross. The scale bar of 5" corresponds to 
67.5 kpc at z=2.765 and 46.3 kpc at z=0.66 (IIo=50kms-i Mpc~^ qo=0). 

Figure 3: Contour map of the VLA 8.44GIIz radio observation of MG1019+0535. The contours are drawn at (3,4,5,8, 16, 32, 64)it 
where a = 0.22mJy/beam. There are no negative contours below the —3a level within this 4"x4" box. The extension seen to the 
south of the core component is not real. 



Figure 4: Keck LRIS spectra of the two optical components of MG1019+0535. Note that whereas component A shows several 
strong emission lines at z=2.765, component B has only one moderately strong line that is close in observed wavelength to the 
IIeIIA1640 line in A. The solid triangles mark the locations of II/3 and [OIII]A5007 under the assumption that the line at 6192A 
is [OII]A3727 at z=0.66; no line emission is detected at these locations. The turn down in the spectrum of B for A >7300A is not 
real and is due to instrumental problems. 
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Table 1: Observed Radio and Optical Flux Densities 



Object 


«1950, 




V 


F. 


Comments 




•^igso 




MHz 


mJy 




Radio Core 


10''16'"56f815(±0f010), 
+5°49'39^'30(±0^'15) 




8439 


3.0±0.4 


VLA 


East Lobe 


+0^'78, -0^'09 




8439 


25.6±1.0 


VLA 


West Lobe 


-0^'54, +0^'22 




8439 


30.5±1.0 


VLA 


Total Source 






365 


925±31 


Texas Catalog 








1400 


454 


White & Becker 1992 








1490 


360±6 


VLA 








4850 


132±19 


Gregory & Condon 1991 








8439 


59.1±1.5 


VLA 


Component A 


10''16'"56f82 
+5°49'39^'0 


4.61 


xlO® 


0.98x10-3 


i?=23.70±0.04 in 3" dia. 
i?c=23.90 (corr. for emission lines) 






1.36 


iXlO® 


6.4x10-3 


K ^ 20 


Component B 


10''16'"56f87 


4.61 


xlO® 


1.79x10-3 


i?=23.05±0.02 in 3" dia. 




+5°49'40^'2 


1.36 


iXlO® 


7.15x10-3 


Jl'=19.89±0.16 


A+B 




4.61 


xlO® 


3.73x10-3 


i?=22.25±0.02 in 7" dia. 






1.36 


iXlO® 


14.9x10-3 


Jl'=19.09±0.24 
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Table 2: Observed Spectral Lines and Their Fluxes 



Object 



'^obs -^ obs 

A 10~^'^erg/s/cm^ 



Wf 



FWHM^fc, 

A A 



ID 



4584.3±0.3 


8.4±0.6 


268 


13.5 


LyaAl216 

ZLya = 1.6(±0 


4664.8±0.8 


2.3±0.6 


66 


10.6 


NVA1240 


5831.9±0.4 


10.4±0.7 


257 


19.6 


CIVA1549 


6174.9±0.5 


8.5±0.6 


170 


21.3 


HeIIAl640 


7179. 6±1. 9 


4.9±1.0 


93 


27.7 


CIII]A1909 


3903: 


3.8: 




11: 


OVIA1034? 


6192. 0±0. 5 


2.4±0.4 


46 


10.4 


OII]A3727? 
or HeIIAl640? 



1) X 10 



43 



erg s" 



25 



Table 3: Comparison of MG1019+0535 Line Ratios with Other Objects 



Spectral 


MG1019 


Average " 


F10214 


TX0211 


NGC 


Model I^ 


Model 11^ 


Line 


+0535A 


Radio Galaxy 


+4724^ 


-122"' 


1068" 


log[7=-1.8 


log[7=-2.8 


Lya 


0.99/ 


9.75 


<0.24 


1.50 


4.75 


19 


21 


NV 


0.27 


0.48 


1.86 


1.33 


1.32 


0.024 


0.0016 


CIV 


1.23 


1.15 


2.10 


1.67 


1.86 


1.13 


0.01 


Hell 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


cin 


0.58 


0.56 


0.57 


0.50 


— 


0.67 


0.16 


"McCarthy (1993) 














''Shields, 


pers. comm. 















<=Elston et al. (1994) 
''van Ojik et al. 1994, 5" aperture 
"Kriss et al. (1992) 

•'^The Lya/HeIIA1640 ratio is 0.69 if we correct for HeIIA1215 contamination using the Hummer and Storey (1987) 
Case B ratio of HeIIA1215/HeIIA1640 for Ne = lQ'^cm.-^ and T=104K. 
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